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FOREWORD 
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Flight Center, Contract NAS8-20082. This task was conducted 
in response to the requirement of Appendix B, Schedule Order 
152, under the direction of Richard Beranek, Unsteady Aero- 
dynamics Branch. 


ACKNOWLEDGEMENTS 

The authors are grateful to Dr. Max F. Platzer of the 
Research Laboratory, Lockheed-Georgia Company, for his 
constant encouragement, guidance and review of the manu- 
script, and to Mr. Robert A. Lott, Lockheed/Huntsville, for 
his constructive comments. 



ii 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC/HREC D162375 


^ SUMMARY 

Oswatitsch and Keune's parabolic method for steady transonic flow is 
applied and extended to bodies of revolution oscillating in a sonic flow field. A 
Laplace transform technique was employed to derive the dipole solution, and 
the Adams-Sears iterative technique was used in the stability derivative cal- 
culation. 

A computer program was developed to perform the stability derivative 
calculation for the slowly oscillating cone and parabolic ogive. Inputs for 
the program are body geometry, thickness ratio, acceleration constant f and 
the pitch-axis location. 

Sample calculations were performed for the parabolic ogive and circular 
cone. These results are compared with those obtained by using other tech- 
niques and the available experimental data for circular cones. For circular 
cones, the numerical studies show that the present results agree quite well 
with the available experimental static stability data and the calculated pitch- 
damping falls between Wehrend and Yanagizawa's experimental results. These 
comparisons show that the present analysis indicates a favorable trend. Finally 
a possible method of improving the present analysis is proposed. 
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NOMENCLATURE 

r+ 2ik, Eq. (2.10) 
body pitch axis location 

freestream sound velocity 
Euler's constant = 0.57721567 

* * # * 

pitch moment coefficient = pitching moment/q^Q (f ) i 

/ ^ 'I* 

normal force coefficient = normal force/ q Q (i ) 

oo 

stability derivatives due to angle of attack 
stability derivatives due to total damping -in -pitch 


in-phase flutter derivatives 
out-of-phase flutter derivatives 

/ * * \ / * 

pressure coefficient =^P -P^j/q^ 
pressure coefficient due to steady axial flow 


in-phase pressure coefficient derivative 


= _ L _ 

cos 0y 9 6 J i 


out-of-phase pressure coefficient derivative = 


cos0 \ ai 


6 —0 


local normal force coefficient slope 

exponential function 
doublet distribution function 
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modified Bessel function of the second kind of zero order 

* * . * 

reduced frequency = U) i /U 
body length 

Laplace inversion transform operator 

sfc jjc 

freestream Mach number = U /C 

oo oo 

"order of" 

Laplace transform parameter 

2 

body cross-sectional area = 2 jtR (x) 
freestream dynamic pressure 
body radius 

time measure in a body-fixed reference frame = t U /i 

00 

freestream speed 

body centerline motion function defined by Eq. (3.5) 

sonic point location 

body-fixed cylindrical coordinates 

Cartesian coordinates 

acceleration constant, Eq. (2.7) 

ratio of specific heats, 1.4 for air 
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body thickness ratio = R /i = R(l) 

second order correction function, Eq. (2.18) 

crossflow out-of-phase perturbation potential 

dipole coordinate 
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ir 

$ 

4 

9 

A. 

9 

fw 

9 

X 

* 


= 3.14159265 

complete time -dependent perturbation potential, Eq. (2.2) 

axial-flow perturbation potential 

pulsating source potential, solution of Eq. (2.10) 

complex perturbation potential, solution of Eq. (2.4) 

dipole solution in Laplace transform plane defined by 
Eq. (2.13) 

reduced complex perturbation potential defined by Eq. (3.2) 
crossflow in-phase perturbation potential 


Subscripts and Superscripts 


( it ( )" t etc. derivatives with respect to the independent variable 


( ) ,( ) ,etc. partial derivatives 

X XX 


( ) 

( ) 

( ) 

( ) 

( ) 

< >! 

<)* 


(!) 

( 2 ) 

(D+(2) 


R 


derivative with respect to time 
first-order quantity 
second -order quantity 

first-order plus second-order quantity 
real part 
imaginary part 

physical quantity with proper dimension 
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Section 1 
INTRODUCTION 

The basic small perturbation equation governing transonic flow over 
slender bodies is well established. Due to its nonlinearity, however, analyt- 
ical (closed form) solutions are obtained only in a few special cases. In 1955, 
Oswatitsch and Keune (Ref. 1) introduced an approximate solution (parabolic 
method) to steady sonic flow over the front part of a half-body of revolution 
by assuming a linear variation of the perturbation velocity in the mean flow 
direction (<f> ~ constant) for linearization of the transonic equation. Maeder 

et al (Refs. 2, 3 and 4) extended this approach to the entire transonic Mach 
number range and applied the theory to full bodies of revolution as well as two- 
dimensional airfoils. Later Maeder and Thommen (Refs. 4 and 5) and Hosakawa 
(Refs. 6 and 7) improved Oswatitsch and Keune' s parabolic method by introducing 
a correction function accounting for the nonlinearity. The parabolic method was 
then improved substantially by Spreiter and Alksne's Local Linearization method 
(Refs. 8 and 9), which yields excellent agreement with experimental results. 
Recently, the method of parametric differentiation by Rubbert (Ref. 10) and 
Rubbert and Landahl (Ref. 11) greatly advances the state of the art of the tran- 
sonic flow problem. Unfortunately, this solution does not appear in a workable 
form and immediate application to the unsteady case is not obvious. 

In the unsteady transonic flow, Landahl (Refs. 12 and 13) following the 

fundamental development of Lin et al (Ref. 14) linearized the unsteady flow 

2/3 2 

equation by assuming high-frequency (k » ( ' for planar flow and k » f In£ 
for axisymmetric flow) oscillation and extensively explored the equations for 
various two and three dimensional shapes. Hsu and Ashley (Ref. 15) further 
studied the blunt-nose slender bodies performing lateral oscillation by extend- 
ing the work of Landahl (Ref. 12). 
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Within the frame of unsteady linearized transonic flow Liu (Ref. 16) 
extended Platzer and Hoffman's (Ref. 17) quasi-slender body theory for 
supersonic flow by adopting a body-fixed coordinate system to treat the half- 
body of revolution oscillating at sufficiently high frequencies in a sonic flow 

field (M = 1). 

00 

In the recent work of Teipel (Ref. 18) and Liu (Ref. 19) they included the 
"nonlinear term" of the unsteady perturbation potential equation in the analysis 
of airfoils (Ref. 18) and bodies of revolution (Ref. 19) oscillating in a sonic flow 
field by extending Oswatitsch and Keune's parabolic method (Ref. 1). Calculated 
results for a NACA 65 A 005 airfoil oscillating at both high and low frequencies 
have been presented by Teipel (Ref. 18). The results showed remarkable agree- 
ment with Nelson and Berman's experimental data (Ref. 20). For the oscillating 
bodies of revolution Liu (Ref. 19) obtained the approximate solutions for both 
high and low frequencies (valid in the range 0^ k < 1) but no numerical results 
were presented. 

The material presented in this report is a continuation of Liu's work 
(Ref. 19). Following his analysis, a brief description of the derivation is 
repeated in this report for the reader's convenience. Stability derivative 
calculations were performed for slowly oscillating parabolic ogives and cir- 
cular cones in a sonic flow field by applying Adams-Sears' (Ref. 21) iterative 
procedure. 

Assuming a constant F (see Eq. 2.7) for a particular body of revolution 
implies that the x- derivative of the perturbation velocity in the main flow 
direction (^ x ) is constant or, at least, changes very slowly. For a parabolic 
half-body of revolution this assumption is nearly correct but is not valid for 


The "nonlinear term" actually is the linear term with variant coefficients 
appearing in the unsteady perturbation potential equation (see Eq. (2.4)). ' 
They are called "nonlinear terms" here to signify that they are the coupling 
terms contributed by the nonlinear term in the steady transonic perturbation 
equation to the unsteady crossflow perturbation equation. 
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cones. In spite of this limitation, calculated results for the circular cone, in 
addition to the parabolic ogive, have been presented in this report (see Section 5) 
to show its qualitative trend. 

Due to the lack of experimental data for the oscillating parabolic ogive, 
no comparison between the calculated and experimental results has been made. 
There are only very few test data available for oscillating cones (Refs. 22,23 
and 24). It is apparent that more experimental work in the field of unsteady 
aerodynamics is urgently needed. 

In the subsequent sections, a body-fixed coordinate system (Fig. 1) is 
adopted and all the quantities are non-dimensional except those superscripted 
with an asterisk (*). 
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Section 2 

APPROXIMATE SONIC EQUATIONS AND SOLUTIONS 


The detailed derivation of an approximate solution to the unsteady "non- 
linear" transonic perturbation potential equation has already been presented 
by Liu in Ref. 19. Here, only a brief description of the derivation is presented. 
It has been shown (e.g., Landahl (Ref. 12)) that the transonic small perturbation 
potential $(x, r,0,t) has to satisfy the following equation 


(1-M 2 )$ + -$ +-y$ aQ +$ - 2M 2 $ , - M 2 

oo xx r r 2 90 rr oo xt oo tt 

r 


= M (y+i) $ $ (2.1) 

oo x xx 


Assuming the body performs simple harmonic motion, the complete perturba- 
tion potential can be written in the form. 


_ a : i_f 

$(x, r,0,t) = <f> (x, r) + <p(x, r, 0) e 


( 2 . 2 ) 


Letting = 1.0 and substituting Eq. (2.2) into Eq. (2.1) yields the following 
equations : 


O(l): * rr .+-* r = ()-+l)* x *„ 


(2.3) 


o<a 0 ): 5> rr +|5> r +'-7? 68 -(2ikJ- x -k 2 ^) = (y+i)^(^ x )_ 


(2.4) 


A A / xkt\ 

where the term < / ? x < l? xx y e j is neglected assuming small amplitude oscilla- 
tions (Ref. 18). 


The basic dipole solution (f can be obtained from the source solution 9? 
by the following relation 

A d<P 

V = 97 cos0 (2.5) 
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By substituting the above relation into Eq. (2.4), one obtains 


+7 <P T - 2lk V x + ^ = O'+D li/a + 


( 2 . 6 ) 


Oswatitsch and Keune (Ref. 1) solved Eq. (2.3) by approximating 


(y+1) 4 > xx = F = constant >0 


(2.7) 


Equation (2.3) becomes 


h +-0 = rt 

rr r r x 


( 2 . 8 ) 


The parabolic solution of Eq. (2.8) has been given by Oswatitsch and Keune 
(Ref. 1) as 

x 

2 -i 


* <x - r) = -57 / f <«) • 


exp 


T] 


4(x-£) 


d£ 

x “£ ‘ 


(2.9) 


Ignoring the term ( y + 1 ) $ x </? in Eq. (2.6) and applying the relation in Eq. (2.7),. 
one reduces Eq. (2.6) to 


1 2 
(p + — (0 - AiO + k (0 = 0 
rr r r r x ■ T 


( 2 . 10 ) 


where A = P + 2ik. 


Laplace transform in the x-direction transforms Eq. (2.10) into 


+7^ - X 2 <^ = 0 , 

rr r r 


( 2 . 11 ) 


For convenience, P is' called the acceleration constant in the subsequent 
sections. 
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where 


^ = (r + 2 ik) p - k^ 


</>(P»r) = 


00 

/ 


(p(x, r) e“P X dx . 


The solution of Eq. (2.11) is 


^(p.r) = -yj F(p) K o (Xr) 


where K Q (Xr) is the modified Bessel function of the second kind of order 


The dipole solution then is 


. 3(P 

^(p,r,0) = g-^cosG = 


^ F(p)X Ko (Xr) 


XF(p) Kj (Xr ) 


Expanding K^(Xr) for small r and performing the inverse Laplace 
form, one obtains the following dipole solution in the physical plane. 


A 


cp(x, r, 0) 



= ZM 

27T r 


COS0 


+ - -§ *** jr (x) [fmy + rj(r, k)J - /V(*) fn(x-0 d^J 

+ - k ^° S9 |f(x) [in^- + V(T, k)J - J** F"(l;).£n(x-t)d^ 

“ ~' Vr S9 j F(x) [ in T~ + k [| “ 2 f F(£)in(x-£)d* 

o 

+ f f"($) (x-|) £n(x - £) d$J, 
o 


6 


(2.12) 

zero. 


. (2.13) 

trans- 


. (2.14) 
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where 


k) = |in(k 2 +^)+c 


1 + i 


‘[§--“' 1 (£)] 


'(2.15) 


Equation (2.14) is the general parabolic dipole solution for Eq. (2.10). Since 
the present report is only concerned with the low frequency case, Eq. (2.14) 
can be further simplified by applying the following "low frequency" assumptions 
(Ref. 19), i.e., 


k < 1 

k = o [rf . 


(2.16) 


Introducing the low frequency parameter JU, = k /P < 1 and requiring the quad- 
ratic and product terms of fl and k be negligible in the expression (2.14) gives 

<P(x, r,0) = + £(x, r;F)]cose (2.17) 

where 

£(x,r ;D = ^If(x) ^n^+c-1^ - J f"(|) in(x - $) d$| 

' o 

+ j F (x) ^£n + cj - J f"($) fn(x - £) d*j . (2.18) 

o 

It has been noted previously (Ref. 19) that the above low frequency solution is 
most suitable for applying the Adams-Sears iterative procedure (Ref. 21). 

Hence, the stability derivative calculation can yield a finite value for the 
limiting case of k — 0. 


*k =0(D=$lim k /r < oo : 0 < k < 1 

r-o 

k = 0 (D 41 im k/r = 0: /JL < 1 

r-o 
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Section 3 

BOUNDARY CONDITIONS AND ITERATIVE PROCEDURES 


This section describes the application of the Adams and Sears iterative 
procedure (Ref. 21) to determine the dipole distribution F(x) in Eq. (2.17). 
The present report follows closely the procedure used in Ref. 16. 

In addition to the boundary conditions at infinity the solution must also 
satisfy the tangency conditions at the body (Ref. 17): 


•= R'(x).(l+^ x ) 
* = -1 + R'(x)*^ 


at r = R(x), 


X r ■ 


Qx-a) + R(x)»R , (xJJ + R (x) • A x ] 


(3.1) 


where 

^(x, r): in-phase perturbation potential, and 
A(x,r): out-of-phase perturbation potential. 

The reduced potential X(x, r) is defined as 

A 

X(x,r) = “5^ = 6^(x,r) + 6X(x,r) . (3.2) 

From the relations in Eqs. (3.2) and (2.17), one can easily obtain the following 
relation 


X(x, r) = |^j- + £(x, r;D . 


(3.3) 
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By differentiating Eq. (3.2) with respect to r and utilizing the relations in 
Eq. (3.1), one gets 


r=R(x) 


= -W(x) + R'(x) • X 


(3.4) 


where 


W(x) is the body center-line motion function and expressed as 
W <*>|r=RM * »+[?-*> .+ §3 «• 


(3.5) 


Applying the Adams -Sear s' iteration procedure one can write 


and 


X(x.r) 

F(x) = F (1) + F (2) 


(3.6) 

(3.7) 


The first term on the right-hand side of Eqs. (3.6) and (3.7) is the slender 
body term and the second term represents the correction for the body thick- 

3|C 

ness and the acceleration constant i . Substituting the above relations into 
Eq. (3.3) and collecting the terms of like order gives: 

.( 1 ) 


,(D _ El 


2ir r 

.( 2 ) 


* (2) = f^ + ^) (x , r;n 


(3.8) 


(3.9) 


From Eqs. (3.4), (3.8) and (3.9) one obtains the following relations: 


2jr r 


_ + = -W (2) (x) + Rf(x). X^ 


at r = R, 


2?r r 


(3.10) 


ft 4 

i is also thickness dependent, but this is related to the nonlinearity of the 
sonic flow. 


9 

LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC/HREC D 162375 


where 


W (1) (x) 


= 6 + (x - a) 6 


W* 2 ’(x) 


2jr 


and 


W (x) = W (1) (x) + W (2) (x) 


(3.11) 


By substituting the relations in Eq. (3.11) into Eq. (3.10) and after a simple 
rearrangement, one can write: 

F (1) (x) = 2 Q(x) fs + (x - a) sj 
- F (*) w = [. 9 ^m] 6 + [. (x . a) SL^l] i + 2 Q (x) !g} 


(3.12) 


(3.13) 


Substituting Eq. (2.18) into Eq. (3.13), F^(x) becomes (see Appendix A for 
details): 


F (2) (x) = 6 


Q ' .fcl +2Z p Q(x) 
It , -K 


+ 6 


(x-a) + 2 Z t Q(x) 

IT 1 


(3.14) 


Upon comparing the relations in Eqs. (3.2) and (3.3) and utilizing the expres- 
sions for F^ 1 ^ (Eq. 3.12), F^ (Eq. 3.14) and ' (Eq. A. 2), one finally 
obtains the in-phase and out-of-phase cross flow potential valid to the second- 
order as: 


*(l) + (2) = S^l £ + Z R ] - ^ + r [r R + Y r ] 


(3.15) 


a ( 1)+(2) = . Qlx) |7 x -a) + Z J - (x-a) + r |v r + yJ , (3.16) 


jrr 


10 
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where 


V R = h °< X > ’ 


I = 4 ^ [ (X " a) Q,(x) + Q(x ll + T? * 
r = £ [(c-DQW-y* jmx-$) <1*1 , 

l—O “ 1 

I = ^ |( c " !) jj x “ a ) Q ( x ) + Q ( X |J jT^ - a ) Q 7 (£) + 2Q'(S) 

+ £ jcQ'(x>- y X Q''($)in(x-$) d$j, • 


(3.17) 


in(x - £) d£ 


Z R V R 


Z I = V I 


2 + in 


r 


4 + Y R ’ and 


2 + in 


/V 


+ Y I • 
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Section 4 

CALCULATION OF STABILITY DERIVATIVES 


The technique of using the results obtained in the preceding sections 
for the calculation of the stability derivatives has been shown in Ref. 17 and 
also used in Ref. 16. Once the in-phase (>&) and the out-of-phase (X) potentials 
are obtained the calculation of the stability derivatives is straightforward. 

4.1 PRESSURE COEFFICIENT 


The pressure coefficient may be expanded for small frequency and am- 
plitude, just as the perturbation potential $ was expanded. That is 


C (x,R,0,t) = 

Jr 


C (x, R) + 


C (x, R) + 6 C (x, R) 

P 1 p 2 -J 


COS0 


(4.1) 


where 



pressure coefficient due to steady axial flow , 
pressure coefficient due to in-phase cross flow, and 
pressure coefficient due to out-of-phase cross flow. 


At = 1, the in-phase and out-of-phase pressure coefficients eval- - 
uated on the body surface r = R(x), valid to the second order, according to 
Platzer and Hoffman (Ref. 17) are 


C = - 2 * (1)+(2) . R' 2 (x) . * (1) 
Pi x x 


= - 2 (* x + *) (1)+(2) + 2 


* (1) + R(x) 




( 1 ) 


- R 2 (x) 


{1) -2R(x) 


(4.2) 
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, and 


ITT 


A^(x, r) = (x-a) 


21*1 

irr 


(4.3) 

(4.4) 


The axial perturbation velocity on the body surface according to Spreiter 
and Alksne (Ref. 9): 


♦(‘Wl = SLM toC jP J s l > r a{ 

x 4 it 4 jr x 4 7 T x - £ 


(4.5) 


The local linearization method by Spreiter and Alksne (Ref. 9) further 
improved the cone-cylinder case by applying the sonic point condition at the 
shoulder. Thus, the improved solution for the cone may be used by letting 
i/2 = 1 in Eq. (64) of Ref. 9: 


<t> (x, R) = f 2 in ex - —■ injf 2 + 

X * ( (y+1) e e 


( 4 . 6 ) 


Expressions of ^ and X are obtained by carrying out partial differentiation 

X X! 

of the appropriate expressions of and A with respect to x: 

*(D = J3l 

X JfR 


*x )+(2> = ar < 1 +z r' + p:( z r) -;^ +e |( v r) toC f !+ ( Y R ) 


x 2 7 T R 


x (1) = -§- + (X-a) 

X fl-R 


_Q_ 

?rR 


. x (l)+(2) _ Q_ 

A x " ,rR 


(x-a) + Z. 


Q 

jrR 


1 + 


( z i). 


(x-a) 


Q'Q" Q ' 2 


ff 2 R 2 jt 2 R 


+ R 


(vA *¥ + (?i) 


(4.7) 
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where and Zj are defined in Eq. (3.17), and 

( v r) x = <V R>' = Ti Q " 

W = < v i>' * £ + 

X 

> 

( Y r) x * < Y R> = £ 

( Y l) = < Y I>' = £ j( c - 1) [q'+ Q"+ (x-a)Q"^j-3I 2 + al 4 - I 


+ 7jH'- i 2 ! 


( z r) x = 

R = 


<Z R* 4ir 


» ro q 

Q (fn^ + c + 1) 


>2 n 
T - r 4_ 


(Z T ) = 


£L 

4tt |L- 


(4.8) 


(x - a) Q"+ 2 Q 7 J^fn-^-^ + c + 1 + (x - a) + Q / 


■31- + al. - I_ 
2 4 5 


a . j[ in m + c + 2 ] Q " + . 


Where I^, 1^ and 1^ are defined as follows: (See Appendix C ) 


■/ 

o 

■/’ 

o 

■f 


Q"{$) in(x-$) d£ 


Q '{£) in(x- £) d£ 


£CT($) £n(x - £ ) d$ 


(4.9) 
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4.2 STABILITY DERIVATIVES 


The four flutter derivatives are defined as (Ref. 17): 



(4.10) 

(4.11) 

(4.12) 

(4.13) 


and related to stability derivatives by (Ref. 25) : 


II 

7? 

U 

Cn 6 

C 

m 

a 

Cm 6 

C N. 4 C N * 
a q 


C ,. + C 
m • m 

a q 

Cm 6 


(4.14) 


In order to be consistent with the slender body assumption, only the 
terms up to second order in 6 are considered in the force and moment coeffi- 
cients. Hence, the last term on the right-hand side of Eqs. (4.11) and (4.13) 
is further simplified in the actual calculation. These two equations may be 
rewritten as . 


15 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



EMSC/HREC D 162375 


™ . _ _ _ . 


" f l 

r 

C ' 
m 6 

IT 

~ Q(l) 

1 (x - a) R Cp dx + i 

' R 2 R C dx 
P3 



o o 


------ 


■ r 1 i 

(•1 

■ :■ . C 

m 6 

7T 

" Q(i) 

y (x - a) RC p dx + I 

R 2 R'C dx 
P 4 

- 


° o 



where 

C = -2 

P 3 x 

C = - 2 (x^> + ^< x >) . 

The expressions of C and C used in this study for a body of revolution 

P 3 P 4 - 

with a radius of R(x) = €Z(x) are given in Appendix B, Eqs. (B.6) and (B.7). 


4.3 CALCULATION PROCEDURE 


To calculate the stability derivatives for a certain pointed body of revo- 
lution, the integration may be carried out according to Eqs. (4.10) to (4.13) 
either analytically or. numerically. In this report the numerical method is 
chosen in order to develop a more general computer program which can handl 
the various body shapes, that are involved. 

To carry out the integration numerically, however, one needs to eval- 
uate the values of the integrand along x coordinate first. To calculate the 

local values of C and C on the body, one first must carry out the integra- 

1 2 

tions that appeared in Eqs. (A.4), (A. 6) and (4.8). These are integrated analyb 
ically as shown in Appendix C. 
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4.4 EVALUATION OF P 

To carry out the calculation of the stability derivatives according to 
those equations obtained in the preceding articles, the value of the accelera- 
tion constant P must be first provided. The methods employed in this report 
are as follows: 


# Circular Cone: 


Based on the assumption Eq. (2.7) and the small perturbation velocity 
for the cone, Eq. (4.6), the following relation can be obtained: 



« 2 (y+i) 



4x 

(y+1) £ 4 e c + 4 (1 -X 2 ) 


(4.17) 


The value of /'evaluated at x = 1/2 for cones of different thickness is plotted 
as curve (T) in Fig. 2. 

• Parabolic Ogive 


Maeder and Thommen (Ref. 4) and Hosokawa (Ref. 6) introduced a method 
to solve for the sonic point on the body and the P value by using the following - 
relations 


* x <x s> = 0 ) 

= nr ) 


In Fig. 2, curve ( 2 ) is given by Maeder and Thommen (Fig. 9 of Ref. 4) for 
the parabolic bodies of revolution, curve ($) is calculated according to con- 
ditions in Eq. (4.18) for the parabolic ogive with the flow potential given by 
Hosokawa (Ref. 6), curve (4) is calculated according to the following relation 
(letting f = 2 in Eq. (50) of Ref. 9): 

P = 2(y +1 ) 7T £ 2 (4.19) 
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Curve ( 2 ) in Fig. 2 given by Maeder and Thommen was originally cal- 
culated for the steady full parabolic body of revolution at Mach one. This 
value of r is used directly for the half parabolic body of revolution in this 
report. Since the physical condition along a parabolic ogive is the same as 
that along the forebody of the parabolic spindle of the same physical thickness 
the sonic point will occur at the same point on both bodies. From Eq. (4.18) 
one can easily see that the same r value can be obtained for both parabolic 
ogive and parabolic spindle of the same thickness. 
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Section 5 

NUMERICAL RESULTS AND DISCUSSIONS 

5.1 APPLICATION TO BASIC BODY CONFIGURATIONS 

To illustrate the effect of body geometry, thickness ratio, pitch axis 
location and the "nonlinear" contribution, i.e., terms associated with T" on 
the static and dynamic stability derivatives, numerical examples were calcu- 
lated for two basic configurations. They are: 

1. Right Circular Cone R(x) = £x 

2. Parabolic Ogive R(x) = cx(2-x) . 

Most of the numerical studies were compared with the previous linear- 
ized flow results (Ref. 16) and other available theoretical and experimental 
results (Refs. 12, 23 and 24). 

5.2 DISCUSSION OF RESULTS 

Stability derivatives were computed and plotted versus three parameters: 
body thickness, pitch axis location and the reduced frequency. In this present 
low frequency calculation, the acceleration constant F has to be preselected 
before the stability derivative calculation. The /"’ values are based essentially 
on three previous steady transonic flow theories, i.e., Maeder-Thommen's 
theory (Ref. 4), Spreiter -Alksne' s theory (Ref. 9) and Hosokawa's theory (Ref. 6). 

In Fig. 2, the acceleration constant P is plotted as a function of thickness ratio 
for the cone and ogive according to the different formulas provided by these theories 

The local normal force slope dC^ /dx for the cone and the parabolic 

a 

ogive are shown in Figs. 3 through 5. In Figs. 4 and 5, the results calculated 
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for the ogive using the value according to Maeder-Thommen and Spreiter- 
Alksne are presented. Since the F values are approximately equal, the results 
obtained using curve ( 3 ) (Hosokawa) are not very different from those obtained 
using curve (2 ) and these results (Hosokawa) are not presented in Figs. 4 and 5. 

In Fig. 6, the normal force coefficients based on different T"’ values 

a 

are plotted against body thickness ratio for a parabolic ogive. Large distinc- 
tion and different trends of the thickness dependence are found for when 

a 

compared with the slender body and the previous linearized flow results. Good 

agreement with experimental results (Ref. 23) for a 12.5-degree cone is found 

in Figs. 7 and 8. In fact, the present calculation predicts, in general, closer 

results to experiment for both CL, and C than the linearized flow results. 

N m 

a a 

This result justifies that the present analysis is a more consistent one. In 

Fig. 9i the damping -in -pitch normal force coefficient C . for a 10-degree 

0 

cone is given in contrast to the various oscillating supersonic cone theories 
(Refs. 26, 27 and 28). 

The damping-in-pitch moment coefficients C . are given in Figs. 10 

m 6 

through 16. Figures 10 through 12 illustrate the pitch-axis location effect 

on for a cone and an ogive. When compared with experimental results 

6 

for a 12.5-degree cone (Fig. 11), the present results fall in between Wehr end's 
(Ref. 23) and Yanagizawa's (interpolated between 10-degree and 15-degree 
cones. Ref. 24) results. Figures 13 and 14 compare the thickness effect for 
a cone and an ogive with various pitch-axis location predicted by the present 
method, linearized flow theory (Ref. 16) and Landahl's theory (Ref. 13). It is 
seen that the slender body theory, in general, predicts a higher cone pitch- 
damping than the other theories. 

Figure 15 presents the frequency dependence of damping-in-pitch moment 

coefficient C for parabolic ogives of different thickness. Previous linearized 
6 
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flow and Landahl' s results are also plotted for comparison. Note the in- 
validity of the linearized and Landahl’ s theory when the reduced frequency 
k approaches zero. In Fig. 16, the influence of body shape on the damping- 
in-pitch moment coefficient is demonstrated for a cone and an ogive of 0.1 
thickness ratio. In contrast to the linearized theory, the present method 
indeed predicts a finite value as the reduced frequency k approaches zero. 
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Section 6 

CONCLUDING REMARKS AND RECOMMENDATIONS 
6.1 CONCLUDING REMARKS 

In the preceding sections, the low frequency potential solution based on 
Oswatitsch-Maeder ' s linearization concept is applied and extended for calcu- 
lating the aerodynamic forces and moments acting upon half bodies of revolu- 
tion performing slow pitching harmonic oscillations. Using this solution, the 
present analysis can provide a finite value for the stability derivatives as the 
reduced frequency approaches zero. This value is of engineering interest 
and cannot be obtained through the previous linearization formulation (Refs. 

13 and 16). 

Following the previous work (Refs. 16 and 17), the Adams-Sears' itera- 
tive procedure was used to obtain the in-phase and out-of-phase perturbation 
potentials, thus facilitating calculation of the stability derivatives. The num- 
erical studies showed that the present results agree quite well with the experi- 
mental value of static stability derivatives for the circular cone (Figs. 7 and 8). 
Comparison with the experimental results (Refs. 23 and 24) for the damping-in- 
pitch moment coefficient for the cone also shows that the present analysis 
indicates a favorable trend. These comparisons seem to substantiate the 
approximation made in the present analysis. However, to obtain the low fre- 
quency solution, it was necessary to approximate the governing equation (2.6), 
by ignoring the term 0 X <P i thus the approach cannot be completely justified 
mathematically. It is, therefore, essential to determine the error incurred 
in the present approximation and to further investigate new results based on 
some other improved solutions for comparison. No definite assessment of 
the technique described in this report can be made based on the currently 
available few test data for circular cones, though it shows a favorable trend. 
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In this study, the technique of selecting the acceleration constant, F , 
has been based on some analytical and semi -analytical methods. If the ex- 
perimental data are available, however, the F value can be calculated from 
either the static or the dynamic test data. The validity of the method de- 
scribed in this report can also be checked by calculating the F value from the 
static data and using this F value to obtain the dynamic data. This calculated 
dynamic data is then compared against the dynamic test data. 

6.2 RECOMMENDATIONS 

Among the unsteady transonic flow studies, Hosokawa's (Refs. 6 and 7) 
nonlinear correction theory seems to offer immediate applicability for the 
present problem. His nonlinear correction procedure can account for the 
transonic nonlinearity, the effect of the transonic normal shock and, in the 
case of unsteady flow, the recovery of the term < ^ X (P XX - By properly incor- 
porating the nonlinear correction procedure to the present low frequency 
solution, one can obtain the "correction value" for the in-phase and out-of- 
phase potential, hence the "correction value" for the stability derivatives. 

This approach is now being pursued, and numerical studies will be performed 
for a parabolic ogive and a parabolic spindle, so that the error incurred due 
to the term < t > x ( P XK can be investigated from the results of parabolic ogive and 
the shock effect can be investigated from the parabolic spindle calculation. 

Due to the lack of experimental data the evaluation of the present method 
is only partially successful. Therefore, it is strongly recommended that 
systematic, both static and dynamic, experiments on the body of revolution in 
the transonic flow range be performed to aid the advancement of the state of 
the art in this particular field of study. 
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Fig. 2 - Acceleration Constant for Bodies of Different Thickness. 
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Fig. 3 - Local Normal Force Slope Variation for Cones at sfc 1. 0. 
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Fig. 4 - Local Normal Force Slope Variation for Parabolic 
Ogive of £ s 0.1 atM s 1.0 

° QO 

30 


LMSC/HREC D1 62375 



o.o 0.2 


0-4 0.6 0.8 • l.o 


Axial Distance (x) 


Fig. 5 - Local Normal Force Slope Variation for Parabolic 

Ogives at M =1.0 
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Fig. 9 - Effect of Mach Number on the Change in Normal Force Coefficient due to 
Pitch Rate for a 10-Degree Cone with Pitch Axis Location at a = 0. 
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Fig. 10 - Damping-in-Pitch Moment Coefficient for a 10-Degree Cone at = 1.0 
vs Pitch Axis Location. 
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Fig. 11 
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Damping -in -Pitch Moment Coefficient for a .12. 5 -Degree Cone at M = 1.0 

00 

vs Pitch Axis Location. 


37 



LMSC/HREC D162375 



0 0-2 0.4 0.6 0.8 1.0 

Body Pitch Axis Location (a) 

Fig. 12 - Damping -in -Pitch Moment Coefficient for Parabolic Ogives at 

M__ = 1.0 vs Pitch Axis Location. 
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Fig. 15 - Damping-in-Pitch Moment Coefficient for Parabolic Ogives at M = 1.0 
» oo 

vs Reduced Frequency with Pitch Axis Locations at a 0.4 and 0.6. 
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Fig. 16 - Damping -in-Pitch Moment Coefficient for a Cone and a Parabolic Ogive 
of the Same Thickness Ratio e = 0.1 at M s = 1.0 versus Reduced 
Frequency with Pitch Axis Location at a = 0.6 
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Appendix A 

SECOND ORDER DOUBLET DISTRIBUTION FUNCTION 


• • 


LOCKHEED -HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC/HREC D162375 


Appendix A 

( 2 ) 

The derivation of F v '(x) is straightforward; however, it involves a 
lengthy step to arrive at the final result, Eq. (3.14). For the purpose of 
keeping record of the result obtained, the step-by-step derivation is pre- 
sented in this appendix. 

For convenience Eqs. (3.12), (3.13) and Eq. (2.18) are rewritten here: 


F (1) (x) = 2Q(x) I 6 + (x-a) 6 


6 + (x - a) 6 


F (2) (x) = 


/ 7 

Q (*) 


7T 


6 + 


f(x.r ;H = 


£r 

Sn 


i k r 
4?r 


r* Cs 

F (x) (fn ^-f- + c 


F (x) { in ^— 7 — + c ] 


- (x - a) Q - 6 + 2 Q(x) 
7 r 

.,) ./■ , 


aj_ 

9r 


( 1 ) 


($) fn(x - £ ) d£ 
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F (£) fn(x - £ ) d£ 


From Eq. (2.18) one may write 


- s 


ikr 


' ,(1) (x) + c — 1^ F- (1) (§) in(x-$) d£ 


4ff 


F ,(1) (x) + c 


■)•/ 


o 
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f" ( 1) (4) fn(x - £ ) d§ 


A- 1 


(3.12) 

(3.13) 

(2.18) 

(A.l) 
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then by making use of the relation in Eq. (3.12) 

£ ( 1 ) (x, r ;T) = ^ jfiQW+sjJx-aJQ^x) + Q(xjjjin^ 


+ c - 1 


- ^ f jfi Q W (x) +6 |(x-a) Q(x) + 2Q'(x) j jfn(x - ojd$ 


o 


+ ^|6Q'(x)+6 |(x-a) Q'(x) +Q(xjj||fn 


Fr‘ 


+ c 


ik 

2?r 


tf 


6 |7x- 


6 Q"(x) +6 (x - a) Q"(x) + 2Q'(x) in(x - $) d£ 


f(x]J |i 


= 6 


ft IKS— 


- lj Q (x) Q"(x) • fn(x - £ ) d£ 


- | (x - a) Q (x) + Q(x)j |in~- + cj - J [If - a) Q"(£ ) 


+ 2tf($)]in(x-$) d$j 


+ 6 


+ cj Q'(x) - y X Q"(x)fn(x-^) d|; 


^ j (x - a) Q'(x) + Q(x)J jjn + c - lj 


A-2 
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US- a) Q (£) + 2Q(£) fn(x-£) d£ 




= 6 i 


Q'(x) in^— + 


£ (c-l)Q'(x) y' X Q / (e)in(x-e)de -Ops] 


+ 6r 


| (x-a) Q'(x) + Q(x) ^^“+ 4 “ ( c - !) (x - a) Q'(x) + Q(x) 


-/ x i -a) Q"(£) + 2Q'(g)]in(x.e) d£ 

A ' 


+ |q'(x) in-^- + c Q'(x) - J Q(i) in(x-£) d£j 


= 6 r 


[ toi r)j£ Q ' W j + £ jlc-DQ W 


•« -r 


Q (|) in(x - £ ) d£. 


+ 6 r 


jr 

4jt 


p-a) Q'(x) +Q(x)] + ^£)j 


+ ~ j(c - 1) [Tx-a) Q'(x) + Q(x)j - f* [pa) □''(£) + 2Q'(£)| in(x-£) d£j 
^ o 




t ^ |cQ(x) - y Q"(!) in(x-!)d|| 


A- 3 
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= 6 r 




V + Y 
R R 


+ 5 r 


[ to^-) v I + y i 


(A. 2) 


where 
V 


_£ 

R ~ 4jr 


Q(x) 


Y = - s — 


r 

4tt 


(c - 1 ) Q 


'<*> -/ 


Q (£) jen(x-|)d^ 


v i = ■ a) Q/(x) + Q(x )] + 


(A. 3) 
(A. 4) 

(A. 5) 


Y I = Tn j (C " [f x ‘ a) Q ( x > + Q(x |] ~J + 2Q'(^i| in( x -g)<Jgj 

v o J 

+ |c Q'(x) -jT Q''<e)in(x-£)d£j . 


Note that the terms containing k^S were ignored for this particular case be- 
cause k is considered to be small and use of the following relation 


6 = 6 q e lkt (6 = ik6) 


has been made to eliminate i . 


ip! 

dr 


|{ni ^r) v R + Y R 


+ 6 r 


[K] 


A -4 
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• 

[7 rA 

. 

r 2 

+ 6 

+Y i 

+ 6 r 

- V, 
r I 


= 6 


V R l 2 + fci T- ] + y r 


: ) 


+ 6 


v i l 2 + to£ 5-J + Y i 


where 


8 Z R + 6 Zj , 


Z R * V R l 2 +tof T-) +Y R 


(A. 7) 


z i = v i l 2 + to ^T~ ) + Y i • 


(A.8) 


Finally, one obtains 


F (2) (x) = 6 


- Q-&) + 2 z 


7 r R 


Q(x) 


+ 6 


s' 2 , 


(x-a)£-&l + 2 Z.Q(x) 
7 r i 


A- 5 
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Appendix B 

FURTHER SIMPLIFICATION OF 
PRESSURE COEFFICIENTS 
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Appendix B 


In order to calculate C and C . in accordance with the assumption 

m r m r 

oo 2 

that the body is slender and terms smaller than order e can be neglected, 

the Cp and in the second term on the right-hand side of Eqs. (4.11) and 
1 2 

(4.13) need to be simplified further. The simplified and , by neglect- 
ing the higher order terms than e, are respectively denoted by and . 

For convenience Eq. (4.2) is repeated here and the order of each term 
in e is noted directly below each term. 


c = - 2 + _ R ,2 (x) . 

Pj X ' ' X 

2 

€ * € 


(B.l) 


C = -2 (A +^r)( 1 ) +(2) +2 

p 2 x 


^ (1) + R(x) 
(«+€)•€ 


.( 1 ) „, 2 , , 

<P X -R (x) 


A {1) - - 2R(x) 


(B.2) 


• (e - c - e ) 


It is obvious thatonly the first term on the right-hand side of Eqs. (B.l) and 
(B.2) needs further consideration. The expressions for and A are: 

X X 


B-l 
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*d)+(2) 


Q Q /2 

= ^' 1+y r>-7^ +v r 


2ir R 


/ ° , 2 . 
e • (e +e ) 


(B.3) 


*<l)+(2) 

* x 


^< 1 + y r i+ ^' y r» 


< V R> 
X Jr R 3 


,o,2. 2 

e • (e + e ) e • e 


(B.4) 


3 3 

e e 


3 3 

6 ; € 


Q_ 

7rR 

(X - a) + Yj 

+ Q. 

ttR 

1 + (Y : ) 

_ 1 x_ 

q'q", q' 2 

- Xf- (x - a) - + (V ) . 

3T R Zir R x 

(B.5) 

€ 4 

(€° + e 2 ) 

e • 

(e ° + e 2 ) 

e 3 • e ° e 3 e 3 



By retaining the terms containing e in Eqs. (B.3), (B.4) and (B.5), one finally 

obtains the expressions for C and C as follows: 

P 3 P 4 




- 2 


Q , . 

irR x ” a 



which can be written as 



(B.6) 

(B.7) 


B-2 
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Appendix C 
INTEGRAL TERMS 
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Appendix C 

To simplify the writing in the text some of the integration terms are 
defined in this appendix as follows: 


o 

■>-/' 

o 

O 

■• = / 
o 


X A 

'(e) in(x-£)d£ 


Q (£) £n(x-£)d£ 


£ Q (£) fn(x - £)d£ 


Q'"(£) fn(x-£)d£ 


£ Q'"(£) fn(x-£)d£ 


(C.l) 


If the radius of the body under consideration is written in the following form 


R(x) = e Z(x) = e(Ax + Bx ) , 


then Eqs. (C.l) become 


j 2 

Ij = It € 

2 A 2 . J 2 + 6 AB • J 3 + 4 B 2 • J 4 


T 2 
h = * c 

2A 2 • Jj + 12 AB • J 2 + 12 B 2 * J 3 


j 2 

I 3 = it e 

2 A 2 • J 2 + 12 AB • J 3 + 12 B 2 • J 4 


T 2 
I 4 = 7T e 

12 AB • J x + 24 B 2 • J 2 


j 2 

I 5 = 7T € 

12 AB • J, + 24 B 2 • J 

— ^ o __ 




> 


> (C.2) 


C-l' 
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where 


x 

j n = J e (n ' u £n( x -e)de 

o 


(n- 1) 

L -(jiT) + tex 

j = 0 


C-2 


n 

, n = 1 , 2, 3, 4 
n 
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